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A  computational  study  of  Mach  5  air  over  a  cylinder  with  a  dielectric  bar¬ 
rier  discharge  actuator  was  performed.  The  actuator  was  pulsed  at  nanosec¬ 
ond  time  scales,  and  it  rapidly  added  thermal  energy  to  the  flow,  creating 
a  shock  wave  that  traveled  away  from  the  pulse  source.  As  the  shock 
wave  traveled  upstream,  it  interacted  with  the  standing  bow-shock,  and 
temporarily  increased  the  bow-shock  standoff  distance.  This  phenomenon 
was  also  observed  experimentally  through  phase-locked  schlieren  photog¬ 
raphy.  This  paper  aims  to  reproduce  flow  phenomena  observed  in  the  ex¬ 
periment  using  high-fidelity  computations  in  order  to  provide  additional 
insight  into  the  shock-shock  interaction,  and  subsequent  effect  on  the  cylin¬ 
der,  through  a  reduced-order  phenomenological  model  of  the  actuator.  A 
three-dimensional  simulation  of  the  experiment  was  able  to  accurately  cap¬ 
ture  the  complex  cylinder /tunnel-sidewall  interaction,  and  to  replicate  the 
changes  in  the  flow  produced  by  the  nanosecond  dielectric  barrier  discharge. 

The  results  show  that  the  device  was  very  effective  at  moving  the  standing 
bow-shock  for  a  minimal  energy  budget. 
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I.  Introduction 


A  recent  experimental  study1  of  a  nanosecond  pulse  dielectric  barrier  discharge  (ns-DBD) 
in  a  Mach  5  flow  demonstrated  the  feasibility  of  a  plasma-based  supersonic  flow  controller. 
In  the  experiment,  a  bow-shock  perturbation  on  a  microsecond  time-scale  was  detected 
in  phase-locked  schlieren  visualization.2  A  compression  wave  was  generated,  due  to  rapid 
localized  heating  from  the  DBD,  which  propagated  upstream  from  the  cylinder  surface  and 
interacted  with  the  standing  bow-shock.  This  interaction  temporarily  increased  the  shock 
standoff  distance,  with  the  series  of  events  repeated  at  an  interval  of  10  p.s  (100  kHz). 

Previous  demonstrations  using  the  ns-DBD  have  included  separated  flow  reattachment  in 
air  flows3  up  to  Mach  0.85,  characterization  of  compression  wave  propagation  in  a  quiescent 
air,2  and  visualization4  of  large-scale,  spanwise  vortex  over  the  airfoil  at  Mach  0.3.  The  flow 
control  mechanism  in  these  experiments2-4  appears  to  be  consistent  with  a  localized  arc  fila¬ 
ment  plasma  actuator  (LAFPA).5-'  This  device  disrupts  the  flow  with  high  amplitude,  high 
bandwidth  perturbations,  with  a  modulation  frequency  near  one  of  the  flow  instability  fre¬ 
quencies,  thereby  triggering  subsequent  growth.  Prior  flow-control  studies,5-7  using  LAFPA 
actuators  in  atmospheric  pressure  jet  flows  for  Mach  0.9  to  Mach  2,  demonstrated  significant 
localized  heating  and  repetitive  shock-wave  formation  by  the  plasma,  large-scale  coherent 
structure  generation,  and  mixing  enhancement.  These  effects  were  achieved  at  a  low  actua¬ 
tor  power  (~10  W  per  actuator),  and  by  forcing  frequencies  near  the  jet-column  instability 
frequency. '  This  low  power  budget  contrasts  with  previous  bow-shock  control  studies,  which 
typically  have  power  budgets  on  the  order  of  10  kW,  using  pulsed  DC  discharges,8  pulsed 
microwave  discharges,9  and  laser  optical  breakdown.10,11 

The  present  work  complements  the  recent  experimental  study1  by  replicating  the  demon¬ 
strated  effect  of  the  ns-DBD  using  the  fluid  dynamics  code,  LeMANS. 12,13  This  was  accom¬ 
plished  by  imitating  the  rapid  energy  coupling  of  the  ns-DBD  via  phenomenological  volumet¬ 
ric  energy  deposition  model.11  The  approach  was  based  on  the  assumption  that  the  primary 
flow  control  mechanism  of  the  ns-DBD  is  rapid  thermal  energy  release.  Many  parameters 
of  interest,  such  as  flow-field  velocity  vectors,  surface  pressure,  and  thermal  loading  were 
not  available  from  the  experiment.  Thus,  the  computational  component  of  the  joint  study 
provided  additional  details  and  insight  about  the  baseline  flow  and  the  flow  perturbation 
mechanism.  In  addition,  the  present  work  developed  a  reduced-order  model  that  replicated 
the  effect  of  the  ns-DBD.  The  reduced-order  model  can  be  used  in  future  studies  to  evaluate 
the  effectiveness  of,  and  optimize  the  use  for,  the  ns-DBD  in  practical  applications,  such  as 
flow  control  in  a  hypersonic  inlet,  isolator,  and/or  engine  exhaust. 
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II.  Experimental  Facilities 


The  experiments  were  conducted  in  a  small-scale  Mach  5  nonequilibrium  wind  tunnel 
that  was  operated  using  dry  air  supplied  from  high-pressure  cylinders,  at  plenum  pressures 
of  p0  =  49.3  kPa,  and  a  mass  flow  rate  of  7  g/s.  The  steady-state  run  time  at  constant 
static  pressure  in  the  supersonic  test  section  was  up  to  10  seconds,  which  was  sufficiently  long 
for  the  blow-down  flow  to  achieve  equilibrium  conditions  within  the  inviscid-core  of  the  test 
section  and  to  conduct  the  ns-DBD  experiments.  The  flow  expanded  through  a  nozzle  with 
a  design  Mach  number  of  5  and  had  a  throat  height  of  1.6  mm,  as  seen  in  Fig.  1.  The  top 
and  bottom  walls  of  the  supersonic  test  section  continue  to  diverge  after  the  nozzle  exit  (at  a 
angle  of  1.5°),  in  order  to  provide  boundary-layer  relief.  Rectangular,  UV-grade,  fused-silica 
windows  were  flush  mounted  on  all  four  walls  in  the  test  section  to  provide  optical  access 
for  schlieren  photography,  emission  spectroscopy,  planar  laser- induced  fluorescence  (PLIF) 
visualization  /  thermometry,15,16  molecular  tagging  velocimetry  (MTV)17, 18  diagnostics,  and 
picosecond  coherent  anti-Stokes  Raman  scattering  (CARS)  measurements.19 


Figure  1.  The  computational  domain  of  a  Mach  5  nonequilibrium  wind  tunnel  with  only  1/4 
of  the  geometry  illustrated  due  to  symmetry. 

A  bow-shock  was  generated  within  the  test  section  by  a  6  mm  outer  diameter  (4  mm  inner 
diameter),  cylinder  model  composed  of  quartz,  which  acted  as  a  dielectric  for  the  ns-DBD 
actuator.  The  model  was  centered  145  mm  downstream  of  the  throat  (35  mm  downstream 
of  the  end  of  nozzle  contour),  where  the  tunnel  cross-section  was  40  mm  x  46  mm  (width  x 
height).  The  model  was  40  mm  long  with  its  ends  affixed  to  the  sidewall  windows. 

Figure  2  provides  a  side-view  schematic  of  the  electrode  configuration  for  the  ns-DBD 
with  the  various  components  proportionally  scaled.  One  actuator  electrode  was  composed 
of  a  10  mm  long,  3  mm  diameter  copper  tube  with  a  wall  thickness  of  0.35  mm,  and  was 
immersed  inside  the  quartz  tube.  Since  the  interior  diameter  of  the  quartz  cylinder  was 
4  mm,  the  tube  electrode  was  positioned  to  contact  the  dielectric  near  the  upstream  side 
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of  the  quartz  cylinder  so  the  plasma-column  induced  during  the  ns-DBD  discharge  was 
approximately  centered  along  the  stagnation  line.  The  other  electrode  consisted  of  a  strip  of 
adhesive  copper  tape  1.5  mm  wide  by  12  mm  long,  and  was  affixed  to  the  windward  surface 
of  the  model.  The  ends  of  the  exposed  copper  tape  electrode  were  covered  by  non-conducting 
polyimide  tape  such  that  the  two  actuator  electrodes  overlapped  over  a  spanwise  distance 
of  10  mm,  with  the  surface  electrode  centered  behind  the  inviscid-core  bow-shock.  Output 
pulse  voltage  and  current  were  measured  during  each  run.  Measured  peak  voltage  and  current 
were  27  kV  and  70  A,  respectively.  Thus,  each  discharge  coupled  between  4-7  mJ/pulse,  and 
lasted  on  the  order  of  100  ns. 


postive  negative 

electrode  electrode 


Figure  2.  Cross-section  of  a  6  mm  quartz  cylinder  with  the  surface  and  immersed  electrodes 
in  the  ns-DBD  experiment. 

Time  evolution  of  the  shock  generated  by  the  ns-DBD  was  captured  on  phase-locked 
schlieren  images,  and  showed  that  the  discharge  pulse  generated  a  compression  wave  which 
propagated  upstream  and  locally  ‘pushed’  the  steady-state  bow-shock  away  from  the  cylin¬ 
der.  The  interaction  temporarily  increased  the  bow-shock  standoff  distance  by  up  to  25%. 1 
The  perturbed  region  bent  away  from  the  flow  stagnation  line,  convected  downstream,  and, 
eventually,  the  shock  envelop  returned  to  its  baseline  shape  about  20  ps  after  the  discharge 
event.  Image  sets  from  the  microsecond-scale  shock-shock  interaction  were  collected  both 
for  a  single-pulse  mode  (pulse  repetition  rate  of  200  Hz),  and  double-pulse  mode.  The  later 
corresponds  with  two  pulses  separated  by  a  10  microsecond  delay  (pulse  repetition  rate  of 
100  kHz).  Additional  details  about  the  Mach  5  ns-DBD  experiment  are  available  in  Ref.  1. 

The  test  section  static  pressure,  =  160  Pa  (±7  Pa),  was  measured  using  a  wall  pressure 
tap  located  on  the  tunnel  sidewall  about  40  mm  upstream  of  the  cylinder  model.  The  pressure 
was  assumed  to  be  constant  through  the  sidewall  boundary,  a  valid  assumption  for  a  thin 
laminar  boundary-layer  (no  experimental  measurements  taken  thus  far  have  indicated  the 
wall  boundary-layer  is  turbulent).  The  freestream  velocity  was  measured  using  nitrogen 
dioxide  (NO2)  MTV.17,18  The  freestream  temperature  was  measured  using  nitric  oxide  (NO) 
PLIF  thermometry  with  a  resultant  temperature  of  T^  =  50  K  to  60  K.  These  measurements 
were  consistent  with  a  freestream  temperature  obtained  using  the  isentropic  flow  relations 


4  of  20 


and  the  freestream  pressure  (T^  =  56  K). 

The  freestream  Mach  number  inferred  from  the  plenum  pressure  and  freestream  static 
pressure  was  Moo  =  4.6,  whereas  a  freestream  Mach  number  of  M^  =  4.8  was  obtained  using 
the  Rayleigh  pitot  formula  and  one-dimensional  normal-shock  relations.  The  freestream 
density  was  inferred  using  the  ideal  gas  relation  (poo  =  p ^  R  Too),  with  the  freestream 
dry  air  composed  of  78%  nitrogen  (N2)  and  22%  oxygen  (02)  by  density.  Table  1  lists  the 
nominal  freestream  conditions  and  uncertainty  bounds  for  the  tunnel’s  test  section. 

Table  1.  Freestream  conditions  and  uncertainty  bounds  for  Mach  5  air  flow  around  a  cylinder. 


Parameter 

Value 

Moo,  [m/s] 

719  ±6 

Too,  [K] 

56  ±5 

Poo,  [g/rn3] 

9.9  ±1.3 

Poo,  [Pa] 

160  ±6.7 

Moo 

4.76  ±0.25 

Figure.  3  shows  a  top-down  view  of  a  schlieren  image  during  a  typical  run  (without  the 
discharge).  The  shock  standoff  distance  was  measured  to  be  AS  =  1.2  mm  and  had  a 
spanwise  length  of  about  10  mm  (25%  of  the  test  section  width).  As  seen  in  the  image,  the 
bow-shock  was  uniform  across  the  span  in  the  inviscid  core  flow. 


Figure  3.  Schlieren  image  of  the  bow-shock  for  Mach  5  air  flow  over  a  6  mm  diameter  cylinder 
(top-down  view). 


III.  Numerical  Methods 

Flow-held  results  were  obtained  using  computational  fluid  dynamics  (CFD)  to  solve  the 
Navier-Stokes  equations.  The  CFD  computations  were  executed  using  the  Michigan  aerother- 
modynamic  Navier-Stokes  (LeMANS)  code,  developed  at  the  University  of  Michigan.20  Le- 
MANS  is  a  general  two-dimensional/axisymmetric/three-dimensional,  parallel,  unstructured 
finite-volume  CFD  code,  and  has  been  used  previously  in  numerous  studies  of  hypersonic 
hows. 12-14  LeMANS  may  be  employed  with  any  of  three  thermodynamic  models:  perfect  gas, 
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equilibrium,  and  nonequilibrium  thermochemistry.  LeMANS  employs  a  two-temperature 
model  to  account  for  thermal  nonequilibrium,  and  a  standard  finite-rate  chemistry  model 
for  nonequilibrium  chemistry.  Thermal  equilibrium  was  assumed  for  the  present  work  due  to 
the  low  freestream  temperature  and  Mach  number.  LeMANS  is  primarily  used  for  steady- 
state  simulations,  but  can  compute  time-accurate  results  with  first-order  temporal  accuracy. 
The  simulations  were  performed  using  second-order  accurate  spatial  discretization,  and  carry 
double-precision  arithmetic  throughout.  Spatial  and  temporal  studies  conducted  in  a  pre¬ 
vious  paper21  indicate  that  the  level  of  grid  resolution  and  time-step  used  achieved  both 
temporally  and  spatially  independent  results. 

The  ns-DBD  controller  used  in  the  experiment  was,  effectively,  a  thermal  actuator,  since  it 
produced  a  flow  control  mechanism  which  appeared  to  be  consisted  with  a  LAFPA.  However, 
unlike  a  LAFPA,  where  the  discharge  occurred  within  a  channel  embedded  in  the  surface,  the 
DBD’s  exposed  electrode  was  flush-mounted,  so  it  also  generated  a  small  body  force  through 
the  induced  ion-wind.  For  the  present  study,  the  momentum  imparted  by  the  ion-wind 
velocity  was  assumed  negligible  in  the  formation  of  the  compression  wave.  This  assumption 
may  be  reevaluated  in  future  work,  but  it  is  unlikely  that  its  inclusion  into  the  model  will 
strongly  influence  the  formation  of  the  subsequent  compression  wave  given  that  the  ion-wind 
would  induce  a  transverse  velocity  (in  contrast  to  the  upstream  traveling  compression  wave). 
In  addition,  recent  high-fidelity  one-dimensional  computations22  of  the  ns-DBD  showed  that 
the  device  created  a  very  small  region  of  relatively-high  ion  density  and  electric  field  required 
for  the  ion-wind  effect,  and  that  the  very-high  electric  field  only  existed  for  a  short  period  of 
time  before  space-charge  shielding  suppressed  the  ion-wind.  As  such,  acceleration  of  neutral 
particles  through  ion-neutral  collisions  was  very  limited,  and  the  majority  of  the  fluid  motion 
was  associated  with  rapid  thermal  energy  transfer  into  the  translational  energy  mode. 

A  phenomenological  model  of  dissipative  heating  was  used  to  represent  the  ns-DBD 
actuator.  The  model  was  implemented  in  the  Navier-Stokes  equations  by  the  addition  of  a 
source  term,  S,  to  the  right  side  of  total  energy  equation,  as  shown  in  Eq.  (1).  Deposition 
of  all  the  energy  into  the  translational-rotational  mode  was  a  strong  assumption,  but  was 
adequate  for  the  purpose  of  this  study  since  it  was  assumed  the  compression  wave  generated 
by  the  DBD  was  due  to  a  rapid  transfer  of  energy  into  the  translation  energy  mode,  as  shown 
in  previous  work  by  Popov.23  The  translational  energy  equation  is: 

dE 

—  +  V-  ((E  +  p)u-T  -u  +  q  +  E(JSp/isp))  =  S  (1) 

where  E  is  the  total  energy  per  volume,  p  is  the  pressure,  u  is  the  mass-averaged  velocity, 
r  is  the  viscous  stress,  and  q  represents  the  total  heat  flux.  The  species  diffusion  flux,  Jsp, 
and  species  enthalpy,  hsp,  represent  energy  transport  for  each  species,  sp.  LeMANS  assumed 
the  fluid  was  continuous,  Newtonian,  and  used  Stokes’  hypothesis  to  determine  the  viscous 
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stresses,  which  were  resolved  using  Wilke’s  semi-empirical  mixing  rule.24 

A  close  examination  of  the  ICCD  images  of  the  plasma  from  the  ns-DBD  pulse  taken  from 
above  the  cylinder  model  in  Fig.  7  in  Ref.  1,  inferred  that  the  highest  concentration  of  plasma 
existed  just  above  the  cylinder  surface  and  wrapped  around  the  cylinder  for  approximately 
30  degrees  on  either  side  of  the  spanwise  stagnation  line  (azimuthal  angle  6  ±  30°).  It  was 
estimated  that  the  plasma  layer  was  0.05  mm  to  0.12  mm  thick,  with  peak  Joule  heating 
assumed  to  occur  about  0.08  mm  from  the  cylinder  surface.  The  image  also  suggested  that 
the  deposition  energy  was  approximately  uniform  within  the  deposition  volume,  so  a  hyper- 
Gaussian  distribution  function  was  used  to  generate  a  near-uniform  deposition  rate  with 
the  deposition  volume,  which  smoothly  and  rapidly  decayed  to  zero  outside  the  region.  The 
hyper-Gaussian  function  was  selected  because  previous  work21  explored  a  standard  Gaussian 
distribution  function  and  found  that  it  lead  to  an  unrealistically  high  peak  temperature  in 
the  flow.  The  deposition  volume  was  wrapped  around  the  cylinder  geometry  by  expressing 
the  phenomenological  ns-DBD  model  in  cylindrical  coordinates: 


S 

r 

9 

z 


Q 

K  ab  c 


exp 


V ((x  ~  Xcy)2  +  (y-  Vcy)2)  -  rc 


tan  1 


y—ycy 

X-XCy 


(z  -  Zc) 


(2) 


where  Q  is  the  total  energy  deposited  and  K  is  a  constant  such  that  J f 'f™  Sr  dr  dOdz  =  Q, 
for  the  selected  values  of  £,  77,  and  £.  For  the  wind  tunnel  geometry  shown  in  Fig.  1,  the 
cylinder  was  centered  at  xcy  =  145  mm,  ycy  =  0,  zcy  =  20  mm  and  had  a  radius  of  3  mm. 
The  center  of  the  phenomenological  energy  deposition  volume  was  0.1  mm  from  the  cylinder, 
so  rc  =  3.1  mm.  The  deposition  volume  represented  the  plasma  glow  observed  in  Fig.  7  in 
Ref.  1  by  setting  a  =  0.04  mm  and  b  =  30°.  The  spanwise  extent  matched  the  overlapped 
electrodes  by  assigning  c  =  5  mm.  Due  to  its  proximity  to  the  cylinder’s  surface,  the  power 
of  the  radial  direction’s  exponential  function  was  £  =  10,  which  ensured  a  nearly-uniform 
deposition  rate  within  the  volume  and  a  sharp  cutoff  outside  of  it.  Some  plasma  may  have 
existed  above  the  polyimide  tape  and  extended  beyond  the  azimuthal  angle  of  30°,  so  the 
radial  and  spanwise  deposition  density  parameters  were  set  slightly  lower:  r)  =  8,  £  =  8. 


IV.  Two-Dimensional  Computations 

The  top-down  schlieren  image  shown  in  Fig.  3  highlighted  the  two-dimensional  structure 
of  the  bow-shock.  A  shock-aligned  structured  mesh  was  developed  for  a  two-dimensional 
computational  domain  around  the  6  mm  diameter  cylinder,  with  the  stagnation  point  located 


7  of  20 


at  x  =  0.003  m.  A  100  x  100  mesh  was  developed  for  the  first  90°  of  the  cylinder  surface 
with  cell  clustering  in  the  radial  direction  at  both  the  cylinder  surface  and  at  the  bow- 
shock,  while  azimuthal  grid  spacing  was  distributed  to  ensure  orthogonality  of  the  cells  at 
the  surface  and  the  shock.  A  simulation  was  carried  out  using  the  inflow  conditions  listed 
in  Table  1  and  assumed  a  no-slip,  isothermal  wall,  =  300  K.  The  wall  temperature  was 
set  to  room  temperature  because  of  the  experiment’s  short  run  time  and  because  a  previous 
study21  explored  other  boundary  conditions  (including  adiabatic,  fully-radiative,  and  partial- 
slip  walls),  and  found  the  wall  boundary  condition  did  not  significantly  contribute  to  the 
bow-shock  standoff  distance.  Temperature  contours  from  the  computation  were  compared 
with  a  side-view  schlicren  image  of  a  bow-shock  in  Fig.  4(a). 


a)  Cylinder  diameter  =  6  mm 


b)  Cylinder  diameter  =  5  mm 


Figure  4.  Contour  lines  of  temperature  for  Mach  5  air  flow  over  a  cylinder.  The  figures 
include  schlieren  images  from  the  experiment  and  Billig’s  empirical  correlation25  for  the  shock 
envelope  (pink  line  with  diamond  symbols). 

As  seen  in  the  figure,  the  experiment’s  bow-shock  standoff  distance,  ASeXp  =  1-2  mm,  was 
20%  smaller  than  the  computational  standoff  distance,  AScfd  =  1-55  mm.  Also  displayed 
in  the  figure  (as  a  pink  line  with  diamond  symbols),  is  an  empirical  relation  estimating 
the  shock  envelope  developed  by  Billig  et  al, 25,26  which  was  in  good  agreement  with  the 
computational  shock  profile.  The  discrepancy  in  the  bow-shock  standoff  location  suggested 
there  were  inappropriate  assumptions  made  in  the  simulation.  An  investigation  exploring 
the  disparity  was  performed  in  previous  work,21  which  found  the  uncertainty  in  freestream 
conditions  accounted  for  a  3%  variation  in  bow-shock  location,  but  no  other  factors,  including 
rarefaction  effects  due  to  partial-slip  walls,  made  an  appreciable  difference  in  the  bow-shock 
standoff  distance.  Finally,  the  two-dimensionality  was  tested  by  computing  flow  for  a  5  mm 
diameter  cylinder  using  the  same  freestream  conditions.  This  scenario  was  also  performed 
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experimentally  (although  the  ns-DBD  study  has  not  been  carried  out  on  the  smaller  diameter 
model).  As  seen  in  Fig.  4(b),  the  5  mm  cylinder  experimental  results  were  in  much  better 
agreement  with  the  computations  and  Billig’s  empirical  relationship.  Thus,  the  mismatch  in 
shock  location  for  the  6  mm  diameter  cylinder  scenario  was  attributed  to  three-dimensional 
effects  found  in  the  experiment  and  a  calculation  of  the  entire  wind  tunnel  was  required. 

V.  Three-Dimensional  Computations 

A  simulation  of  the  entire  tunnel  (including  the  nozzle  and  the  region  downstream  of  the 
test  section),  was  performed  to  identify  the  reason  for  the  mismatch  in  bow-shock  location 
observed  in  the  two-dimensional  simulations.  Details  about  the  tunnel  (including  the  throat 
dimensions)  were  provided  in  Section  II,  with  Fig.  1  outlining  of  the  computational  domain 
used  in  the  simulation.  Only  one  quarter  of  the  tunnel  was  simulated  because  the  flow  was 
laminar  and  the  tunnel  was  symmetric  in  both  the  spanwise  and  transverse  directions.  The 
nozzle  throat  conditions  (denoted  with  a  superscript  *),  are  listed  in  Table  2. 

Table  2.  Nozzle  throat  conditions  for  the  Mach  5  wind  tunnel. 


Parameter 

Value 

u* ,  [m/s] 

318 

T*,  [K] 

250 

P*,  [kg/m3] 

0.361 

p*,  [kPa] 

26 

M* 

1.0 

A  computational  grid  was  developed  for  the  tunnel  using  four  computational  blocks 
composed  of  structured  cells.  The  grid  spacing  was  such  that  grid  clustering  occurred  near 
all  walls,  with  at  least  20  points  defining  each  boundary  layer.  This  level  of  resolution  was 
sufficient  for  resolving  the  laminar,  wall-bounded  flow,  since  A <  100  for  all  surfaces, 
where  s  represents  the  wall-normal  grid  spacing  along  each  respective  wall,  scaled  by  local 
inner  coordinates  (i.e.,  the  friction  velocity  divided  by  the  kinematic  viscosity).  A  grid 
resolution  study  was  performed  in  previous  work21  for  the  computational  block  surrounding 
the  cylinder  geometry,  which  also  accounted  for  the  unsteady  shock  interaction  due  to  the 
ns-DBD  pulse.  The  grid  surrounding  the  cylinder  was  uniformly-spaced  in  both  the  x  and 
6  directions  to  capture  the  propagation  of  the  ns-DBD’s  compression  wave  and  subsequent 
bow-shock  perturbation.  In  total,  the  three-dimensional  computational  domain  contained 
about  15  million  cells  and  was  run  using  512  processors  (with  a  parallel  efficiency  of  ~75%). 
The  grid  used  in  the  simulation  had  more  than  twice  the  number  of  cells  than  a  previous 
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simulation  of  the  tunnel  performed  by  Nishihara  et  a/.16 

No  slip,  isothermal  walls  (T^,  =  300  K)  were  assumed  for  all  surfaces  and  a  non- 
reflective  first-order  extrapolation  was  used  at  the  domain  exit  plane.  For  the  baseline 
simulation  (i.e.,  without  the  discharge  event),  implicit  time  integration  was  employed  with 
a  time  step-size  varying  from  At  =  0.1  ns  to  At  =  2  p.s.  The  simulation  was  started  from 
quiescent  air,  with  the  inflow  boundary  conditions  specified  in  Table  2.  As  a  result  of  the 
quiescent  starting  conditions,  the  simulation  required  ~10  ms  of  simulation  time  (~20,000 
computational  time-steps),  for  the  wall  boundary- layer  and  cylinder  bow-shock  structure  to 
fully  develop.  The  simulation  was  allowed  to  evolve  for  an  additional  30  ms  of  simulation  time 
(30,000  computational  time-steps  at  At  =  2  p.s),  to  capture  the  low  frequency  oscillation  of 
the  sidewall  boundary-layer/cylinder  juncture. 

Figure  5  highlights  the  boundary- layer  and  bow-shock  structure  for  the  fully-developed, 
quasi-steady,  solution.  As  seen  in  the  figure,  substantial  growth  of  the  sidewall  boundary- 
layer  was  observed  upstream  of  the  cylinder  model,  in  contrast  to  moderate  growth  of  the  top- 
wall  boundary-layer.  The  difference  was  due  in  part  to  the  1.5°  divergence  of  the  top  tunnel 
wall.  Figure  5  also  shows  the  midspan  planar  contour  through  the  top-wall  boundary-layer 
which  shows  the  boundary-layer  temperature  was  appreciably  higher  than  in  the  inviscid-core 
flow  due  to  the  isothermal  walls  and  acceleration  of  the  core  flow. 


Figure  5.  Mach  1.25  iso-contour  colored  by  temperature  at  t  =  19  ms  for  Mach  5  air  flow  in 
the  tunnel  over  a  6  mm  diameter  cylinder  centered  145  mm  downstream  of  the  tunnel  throat. 

Figure  6  shows  schlieren  images  from  the  experiment  which  have  been  overlaid  with 
the  computational  solutions.  The  schlieren  images  were  taken  with  a  knife-edge  set  in  the 
streamwise  direction.  This  was  replicated  from  the  computational  results  by  plotting  con¬ 
tours  of  the  streamwise  density  gradient  (with  gradients  less  than  100  kg/rn4  excluded  for 
clarity).  The  thickness  and  location  of  the  polyimide  tape  placed  on  the  cylinder  was  also 
identified.  As  seen  in  the  figure,  the  computational  bow-shock  location  and  shock  shape  are 
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in  much  better  agreement  with  the  experiment,  especially  compared  to  the  two-dimensional 
results  presented  in  Section  IV. 


Figure  6.  Schlieren  images  of  Mach  5  air  flow  over  a  6  mm  diameter  cylinder  with  the  top 
half  of  each  image  overlaid  by  the  computational  streamwise  density  gradient. 


Aside  from  schlieren  visualization,  the  only  other  experimental  datum  available  was  a 
surface  pressure  measurement  located  on  sidewall  of  the  tunnel,  40  mm  upstream  of  the 
cylinder.  The  computed  pressure  at  the  location  of  the  pressure  tap  was  pw  =  147  Pa, 
which  was  slightly  lower  than  the  experimental  value  of  pw  =  160  Pa.  Given  the  precision 
of  the  experimental  measurement  (±7  Pa),  and  the  uncertainties  associated  with  the  tun¬ 
nel  conditions,  the  computational  results  were  consider  to  be  in  good  agreement  with  the 
experiment. 

In  an  ideal  wind  tunnel,  the  inviscid  flow  monotonically  accelerates  to  the  desired  condi¬ 
tions  in  the  test  section.  However,  the  small  Mach  5  tunnel  used  in  this  work  had  relatively- 
thick  tunnel  wall  boundary-layers.  The  thick  sidewall  boundary-layer,  highlighted  in  Fig. 
5,  existed  far  upstream  in  the  test  section.  This  obstacle  caused  the  inviscid,  supersonic, 
core-flow  area  to  shrink  and,  subsequently,  the  velocity  slowed  upstream  of  the  test  section 
(x  ~  40  mm  upstream  of  the  cylinder).  The  decreased  Mach  number  coincided  with  an 
increase  in  static  pressure  within  the  inviscid  core  and,  a  higher  than  anticipated  post-shock 
pressure.  Figure  7  shows  a  slice  of  pressure  contours  along  the  transverse  symmetry  plane. 
As  seen  in  the  figure,  there  was  a  large  pressure  drop  in  the  spanwise  direction  just  after  the 
bow-shock.  The  strong  pressure  gradient  was  due  to  the  low  pressure  in  the  boundary-layer, 
and  was  exacerbated  by  the  presence  of  a  necklace  vortex  that  formed  within  the  sidewall 
boundary-layer  just  upstream  of  the  cylinder. 

The  large  pressure  gradient  caused  a  significant  portion  of  the  post-shock  inviscid  core 
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Figure  7.  Top-down  view  of  velocity  streamlines  and  pressure  contours  for  the  transverse 
symmetry  plane  ( y  =  0),  for  air  flow  in  a  Mach  5  wind  tunnel  over  a  6  mm  diameter  cylinder. 


flow  to  turn  into  the  spanwise  direction.  While  the  post-shock  flow  escaping  into  the  bound¬ 
ary  layer  was  subsonic,  it  had  a  high  velocity,  so  a  large  portion  of  the  post-shock  mass-flow 
was  not  two-dimensional.  This  three-dimensional  effect  caused  the  bow-shock  location  to 
move  much  closer  to  the  cylinder  than  would  be  expected  in  a  two-dimensional  flow. 

The  formation  of  the  necklace  vortex  was  the  result  of  a  complex  interaction  between 
the  cylinder  and  the  tunnel  sidewall.  The  vortex  experienced  a  slight  oscillation  at  a  fairly 
low  frequency  (30  Hz).  The  low  frequency  oscillation  “ breathing ,”  of  the  boundary  layer  is  a 
characteristic  of  wall/blunt-body  juncture  flows,  where  a  A-shock  structure  forms  at  the  edge 
of  the  inviscid  region  and  the  boundary  layer.2'  The  oscillation  occurred  in  the  streamwise 
direction  and  produced  periodic  changes  to  the  cylinder /sidewall  wake,  which  can  be  seen 
in  Fig.  8.  Oscillations  at  the  shock/boundary  layer  junction  were  also  observed  in  schlieren 
images  from  the  experiment.  In  addition,  Fig.  8  shows  the  existence  of  a  reflected  wave  that 
appears  downstream  of  the  cylinder  wake  along  the  flow  spanwise  centerline,  which  impinges 
on  the  recirculation  zone  in  the  cylinder’s  wake. 

VI.  Energy  Deposition 

With  the  three-dimensional  simulation  results  able  to  replicate  the  experiment’s  baseline 
flow,  and  able  to  provide  evidence  explaining  the  bow-shock  location  mismatch  observed  in 
the  two-dimensional  computations,  a  simulation  of  the  ns-DBD  discharge  event  was  carried 
out.  The  simulation  assumed  a  total  energy  deposition  Q  =  2  kW  for  50  ns  (for  1/4  the 
geometry),  at  a  time  step  of  At  =  0.5  ns  to  ensure  temporal  independence  based  on  previous 
two-dimensional  simulations.21  Thus,  the  total  energy  deposited  was  0.4  mJ/pulse,  which 
corresponds  to  a  rapid  thermalization  of  10%  of  the  input  energy.  This  was  qualitatively 
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Figure  8.  Mach  1.25  iso-surface  colored  by  temperature  for  Mach  5  air  in  the  tunnel  over  a 
6  mm  diameter  cylinder,  illustrating  the  breathing  of  the  tunnel  sidewall/cylinder  wake  over 
a  30  ms  cycle  (~30  Hz). 


consistent  with  kinetic  modeling  calculations,28  which  predicted  energy  thermalization  to 
occur  predominantly  within  a  few  tens  of  nanoseconds  after  the  discharge  pulse,  with  up  to 
~20%  over  the  first  50  ns.  While  Popov23  suggests  the  overall  thermal  efficiency  of  these 
types  of  actuators  may  be  as  high  as  30%,  an  efficiency  of  10%  is  a  reasonable  estimate  for 
the  first  50  ns  of  the  discharge. 

The  simulation  was  run  from  the  baseline  solution  for  17  ps  using  a  At  =  0.5  ns  time- 
step  (i.e. ,  34,000  iterations)  to  capture  the  evolution  of  the  compression  wave/bow-shock 
interaction.  Figure  9  shows  phase-locked  schlicren  images  at  various  times  after  the  ns- 
DBD  pulse  along  with  the  computational  streamwise  density  gradient.  As  seen  in  the  figure, 
the  computations  are  in  excellent  agreement  with  the  experiment  for  the  first  6  ps,  which 
corresponds  to  the  maximum  extent  of  the  experiment’s  perturbed  bow-shock. 

Figure  10  plots  distributions  of  the  nondimensional  pressure  coefficient,  Cp  =  2  (p  — 
Poo)/(Poo  ulo)i  and  nondimensional  heat  transfer  coefficient,  Ch  =  2  q/(poo  along  the 
surface  of  the  cylinder  at  various  times  after  the  ns-DBD  pulse.  The  curves  were  spanwise- 
averaged  over  the  region  of  the  cylinder  subjected  to  the  inviscid  core-flow  (0.015  m  <  z  < 
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Figure  9.  Side-view  schlieren  images  for  Mach  5  air  flow  over  a  6  mm  diameter  cylinder 
with  the  top  half  of  each  image  overlaid  by  the  computational  streamwise  density  gradient  for 
various  time  delays  after  the  ns-DBD  pulse. 


0.025  m),  since  that  region  corresponds  with  the  spanwise  extent  of  the  ns-DBD  deposition 
model  and  was  consistent  with  the  placement  of  the  electrodes  in  the  experiment. 


(b)  Nondimensional  heat  transfer  coefficient. 


Figure  10.  Spanwise-averaged  distributions  of  nondimensional  pressure  and  heat  transfer 
coefficients  along  the  surface  of  the  cylinder  in  the  inviscid  core  (0.015  m  <  z  <  0.025  m). 


As  seen  in  Fig.  9,  the  rapid  energy  deposition  of  the  ns-DBD  creates  a  strong  compression 
wave  that  moves  upstream  to  interact  with  the  standing  bow-shock.  However,  the  rapid 
volumetric  deposition,  which  was  centered  a  small  distance  away  from  the  cylinder,  also 
created  a  streamwise  moving  wave  which  interacted  with  the  cylinder  and,  momentarily,  but 
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significantly,  increased  the  stagnation  region  surface  pressure  (6  <  45°),  as  seen  in  Fig.  10(a). 
The  elevated  surface  pressure  corresponds  with  an  increase  in  drag  on  the  cylinder,  but  the 
drag  augmentation  starts  to  subside  as  the  wave  expands  around  the  cylinder.  About  4  ps 
after  the  pulse,  the  surface  pressure  fell  below  the  baseline  solution,  which  coincided  with  a 
slight  reduction  in  total  drag  on  the  cylinder  versus  the  baseline  scenario.  Just  prior  to  this 
event,  the  upstream- moving  compression  wave  began  interacting  with  the  bow-shock  (~3  ps 
after  the  pulse),  which  caused  the  bow-shock  to  move  away  from  the  cylinder.  As  a  result, 
the  surface  pressure  continued  to  drop  on  the  windward  side  of  the  cylinder  until  ~10  ps 
after  the  pulse,  which  corresponded  with  the  maximum  outward  extent  of  the  bow-shock  in 
the  simulation.  As  the  solution  continued  to  evolve,  the  surface  pressure  gradually  returned 
to  its  original  distribution  and  the  bow-shock  moved  back  to  its  original  standoff  distance. 

The  nondimensional  heat  transfer  distributions  in  Fig.  10(b)  show  similar  behavior,  with 
the  rapid  energy  deposition  process  raising  the  local  translational  temperature  and  causing 
a  large  increase  in  surface  heating.  Surface  profiles  shortly  after  the  ns-DBD  pulse  show  that 
the  peak  heating  did  not  occur  along  the  stagnation  line,  but  rather  was  located  about  20° 
away  from  the  symmetry  plane.  This  was  due  to  the  uniformity  of  the  energy  deposition 
in  the  azimuthal  direction  and  that  the  deposition  was  centered  a  small  distance  from  the 
surface,  so  the  heated  flow  moved  away  from  the  stagnation  line  as  it  moved  toward  the 
cylinder.  Unfortunately,  high  heat  loading  on  the  windward  side  of  the  cylinder  persisted, 
even  after  the  initial  wave  moved  past  the  cylinder  and  the  bow-shock/compression  wave 
interaction  reached  its  furthest  extent.  This  occurred  because  the  energy  was  deposited  in 
the  stagnation  region  of  the  flow,  which  required  a  longer  time  to  dissipate  than  if  it  had 
been  placed  in  the  freestream.  In  addition,  radiation  in  the  visible  spectrum  was  observed 
in  the  experiment,1  but  the  fluid  code  lacked  a  radiation  model  to  remove  energy  from  the 
computational  domain.  Likewise,  the  imposed  isothermal  wall  boundary  condition  did  not 
allow  for  any  of  the  deposition  energy  to  be  absorbed  by  the  cylinder. 

Figure  10  also  shows  that  aft  portion  of  the  cylinder  was  largely  unaffected  by  the  ns- 
DBD  pulse  do  to  the  large,  very-low  density,  recirculation  zone  established  in  the  cylinder’s 
wake.  The  extent  of  the  recirculation  region  was  apparent  in  Fig.  8,  which  shows  the  region 
stretching  about  0.24  m  (4  cylinder  diameters)  downstream  of  the  model.  Overall,  the 
ns-DBD  provided  a  slight  decrease  in  total  drag  over  the  entire  ns-DBD  cycle,  but  the 
improvement  was  dwarfed  by  the  large  increase  in  peak  and  total  heating  experienced  by  the 
cylinder. 

As  previously  noted,  the  computations  diverge  slightly  from  the  experiment  ~6  ps  after 
the  ns-DBD  pulse.  While  the  computational  compression  wave  speed  was  similar  to  that 
observed  in  the  experiment,  its  strength  was  greater.  Thus,  the  perturbed  bow-shock  con¬ 
tinued  outward  for  an  additional  3  ps  and,  consequently,  the  perturbed  bow-shock  also  took 
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longer  to  return  to  its  baseline  state.  These  results  indicate  that  in  the  model,  S,  the  total 
coupled  energy,  Q,  and/or  the  duration  of  the  deposition,  r,  may  need  adjustment.  In  par¬ 
ticular,  the  compression  wave  appeared  overly  strong  near  the  stagnation  line,  which  implies 
that  the  azimuthal  extent  of  the  deposition  volume  should  be  increased  (i.e. ,  b  >  30°).  This 
hypothesis  was  based  on  a  close  investigation  of  Fig.  9  (e),  which  showed  a  slight  kink  in  the 
simulation’s  perturbed  bow-shock  6  ~40°  from  the  stagnation  line.  This  kink  was  also  ob¬ 
served  and  highlighted  on  the  schlieren  image,  though  it  occurred  ~60°  from  the  stagnation 
line.  This  inconsistency  suggests  the  model  volume  did  not  extend  far  enough  around  the 
cylinder  and/or  the  parameter,  77,  which  controlled  the  deposition  decay  rate  in  the  radial 
direction  was  too  large  and  should  be  reduced. 

Expanding  the  size  the  deposition  region  will  weaken  the  compression  wave  strength, 
while  not  significantly  slowing  its  speed.  A  weaker  wave  should  allow  the  perturbed  bow- 
shock  to  return  to  its  nominal  location  at  a  rate  closer  to  that  observed  in  the  experiment, 
since  the  temperatures  in  the  stagnation  region  would  be  lower  and  the  perturbed  bow-shock 
would  not  be  displaced  as  far.  Likewise,  performing  an  additional  simulation  of  the  entire  in¬ 
teraction  using  an  adiabatic  wall  boundary  condition,  or  better,  coupling  the  cylinder  surface 
boundary  to  a  three-dimensional  conjugate  heat-transfer  model,  like  the  model  developed 
by  Anion,29  would  allow  some  of  the  deposition  energy  to  be  absorbed  into  the  cylinder 
surface  as  the  solution  evolved.  This  change  in  the  fluid  code’s  run  conditions  should  also 
allow  the  perturbed  bow-shock  to  relax  back  to  its  baseline  state  more  quickly  since  it  would 
provide  an  additional  avenue  for  the  deposition  energy  located  in  the  stagnation  region  dis¬ 
sipate  out  of  the  solution  domain.  However,  the  extra  computations  needed  to  fine-tune  the 
three-dimensional  phenomenological  model  and  cylinder  boundary  condition  are  beyond  the 
current  resources  dedicated  for  this  work. 

The  results  demonstrate  the  need  for  developing  a  non-empirical,  physics-based  kinetic 
model  of  surface  nanosecond  pulse  dielectric  barrier  discharge  to  reduce  the  number  of  ad¬ 
justable  parameters  in  the  phonological  model.  The  physics-based  kinetic  model  needs  to 
predict  spatial  and  time  distributions  of  the  energy  coupled  and  thermalized  in  the  ns-DBD 
actuator,  in  terms  of  actuator  geometry,  voltage  pulse  shape,  temperature,  and  pressure. 
Such  a  model  has  been  developed  in  Ref.  30,  but  without  air  plasma  kinetics  incorporated. 
Given  the  high  computational  cost  associated  with  evaluating  the  three-dimensional  simu¬ 
lations,  it  may  prove  more  feasible  to  experimentally  explore  the  ns-DBD  on  the  smaller 
5  mm  diameter  cylinder  since  three-dimensional  effects  were  less  pronounced,  or  to  conduct 
additional  repetitions  of  the  existing  ns-DBD  system  for  conditions  with  thinner  sidewall 
boundary-layers,  so  two-dimensional  simulations  of  the  phenomenon  can  be  evaluated. 

While  the  three-dimensional  computation  results  did  not  perfectly  replicate  the  shock 
behavior  observed  in  the  experiment,  they  did  match  the  experiment  quite  well  for  the 
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first  6/is  of  the  ns-DBD  pulse.  In  addition,  the  solution  showed  that  the  compression  wave’s 
interaction  with  the  bow-shock  was  directly  responsible  for  the  bow-shock  movement,  and  not 
the  result  of  a  spanwise  portion  of  the  compression  wave  interacting  the  with  necklace  vortex 
located  inside  the  sidewall  boundary-layer.  Had  the  latter  occurred,  the  interaction  between 
the  compression  wave  and  the  necklace  vortex  would  have  generated  spanwise  relief  to  the 
flow  (i.e. ,  the  flow  would  become  more  two-dimensional).  As  seen  in  the  two-dimensional 
analysis  of  the  5  mm  diameter  cylinder  flow  and  the  discussion  in  the  analysis  of  the  three- 
dimensional  flow,  a  more  two-dimensional  flow  would  also  experience  an  increase  in  the 
bow-shock  standoff  distance  since  less  of  the  post-shock  flow  would  escape  into  the  sidewall 
boundary-layer.  This  was  not  the  case,  since  the  compression  wave’s  interaction  with  the 
core-flow  bow-shock  began  3  ps  after  the  start  of  the  ns-DBD  discharge,  which  was  much 
sooner  than  the  compression  wave’s  interaction  with  the  necklace  vortex.  In  addition,  the 
necklace’s  vortex  unsteady  influence  on  the  flow  occurred  on  a  time-scale  three  orders  of 
magnitude  slower  than  those  observed  using  the  ns-DBD  actuator.  As  such,  the  complex 
sidewall  flow  did  not  have  sufficient  time  to  influence  the  core-flow  perturbations  studied 
here. 


VII.  Conclusions 

A  computational  study  of  Mach  5  air  flow  around  a  6  mm  cylinder  with  a  nanosecond 
DBD  discharge  was  carried  out  using  high-fidelity  numerical  simulations  to  explore  the  im¬ 
pact  the  discharge  had  on  the  flow.  The  shock  standoff  distance  for  the  6  mm  cylinder 
scenario  was  found  to  be  20%  smaller  in  the  experiments  than  the  predictions  of  the  empiri¬ 
cal  correlation  and  two-dimensional  computations.  However,  the  shock  standoff  distance  for 
a  5  mm  cylinder  computation  agreed  well  with  the  experiment.  A  three-dimensional  simula¬ 
tion  of  the  entire  wind  tunnel  was  subsequently  performed,  which  replicated  the  bow-shock 
structure  seen  in  schlicren  photographs,  predicted  the  width  of  the  tunnel’s  inviscid  core,  and 
matched  the  measured  sidewall  static  pressure.  Having  matched  the  experiment’s  bow-shock 
location  in  the  three-dimensional  computation,  the  discrepancy  in  the  bow-shock  location 
observed  in  the  two-dimensional  simulations  was  attributed  to  the  formation  of  a  A-shock 
structure  near  the  junction  of  the  cylinder  and  sidewall.  As  a  result,  the  flow  developed 
a  thick  sidewall  boundary- layer  upstream  of  the  test  section,  and  a  strong  necklace  vortex 
that  formed  around  the  cylinder  at  the  sidewall  junction.  These  features  resulted  in  a  com¬ 
plex,  quasi-steady,  boundary-layer  which  increased  the  post-shock  spanwise  velocity,  thereby 
drawing  the  bow-shock  closer  to  the  cylinder,  even  though  the  experiment’s  schlicren  im¬ 
ages  and  the  simulation’s  streamwise  density  gradient  contours  both  showed  the  bow-shock 
structure  to  have  very  little  spanwise  variation  within  the  inviscid-core  flow  test-section. 
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A  three-dimensional  simulation  of  the  entire  tunnel  with  the  nanosecond  pulse  was  also 
carried  out  and  was  able  to  match  the  compression  wave  speed  and  the  resultant  perturbed 
shock  wave  speed  for  the  first  6  ps  after  the  deposition,  which  corresponded  to  the  peak 
location  of  the  perturbed  bow-shock  in  the  experiment.  However,  the  compression  wave 
in  the  computation  was  too  strong  near  the  stagnation  line  which  pushed  the  bow-shock 
outward  an  additional  20%  over  the  next  3  ps.  While  these  results  do  not  agree  perfectly 
with  the  experiment  during  the  later  portion  of  the  discharge  cycle,  the  behavior  of  the  shock 
perturbation  was  found  to  be  independent  of  the  form  of  energy  input  (due  to  the  short 
pulse  duration),  so  the  resultant  flow  perturbation  was  accurately  simulated.  In  addition, 
the  results  indicate  that  the  cylinder /sidewall  junction  was  not  immediately  influenced  by 
the  discharge  event,  so  the  change  in  the  bow-show  standoff  distance  was  a  direct  result  of 
the  nanosecond  pulse  surface  dielectric  barrier  discharge,  and  not  a  secondary  effect  due  to 
the  compression  wave’s  interaction  with  the  necklace  vortex  in  the  sidewall  boundary-layer. 

The  nanosecond  DBD  actuator  produced  a  strong  compression  wave,  which  yielded  a 
minor  decrease  in  total  drag  to  the  cylinder  over  the  lifetime  of  the  discharge  cycle  at  the 
expense  of  a  large  increase  in  peak  surface  heating.  While  the  device  was  not  very  effective  at 
improving  the  cylinder’s  surface  conditions,  it  was  very  effective  at  moving  a  strong  standing 
shock  for  a  very  minimal  energy  budget.  As  such,  this  technology  could  prove  very  useful 
in  supersonic  inlets  and  isolators  found  in  RAM  jet  and  SCRAM  jet  engines,  where  engine 
‘un-start’  continues  to  be  a  technical  challenge. 
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